THE STRESS INTENSITY FACTOR FOR NON-SMOOTH 
FRACTURES IN ANTIPLANE ELASTICITY. 
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Abstract. Motivated by some questions arising in the study of quasistatic 
growth in brittle fracture, we investigate the asymptotic behavior of the en- 
ergy of the solution u of a Neumann problem near a crack in dimension 2. We 
consider non smooth cracks K that are merely closed and connected. At any 
point of density 1/2 in K, we show that the blow-up limit of u is the usual 
"cracktip" function y/rsm(0/2), with a well-defined coefficient (the "stress in- 
tensity factor" or SIF). The method relies on Bonnet's monotonicity formula 
[2] together with T-convergence techniques. 



According to Griffith's theory, the propagation of a brittle fracture in an elas- 
tic body is governed by the competition between the energy spent to produce a 
crack, proportional to its length, and the corresponding release of bulk energy. An 
energetic formulation of this idea is the core of variational models for crack propa- 
gation, which were introduced by Francfort and Marigo in [10] and are based on a 
Mumford-Shah-type [16] functional. 

In this work, we will restrict ourselves to the case of anti-plane shear, where 
the domain is a cylinder f2 x R, with f2 C R 2 , which is linearly elastic with Lame 
coefficients A and fi. Moreover we assume the crack to be vertically invariant, while 
the displacement is vertical only. Under those assumptions, the problem reduces to 
a purely 2D, scalar problem. Extending our result to (truely 2D) planar elasticity 
requires a finer knowledge of monotonicity formulas for the bilaplacian and is still 
out of reach, it is the subject of future study. 

Given a loading process g : t h-> g(t) £ i/ 1 (fi), and assuming that K(t) C fi (a 
closed set) is the fracture at time t, the bulk energy at the time tg is given by 



where the minimum is taken among all functions u G 1 (17 \ K(to),M) satisfying 
u = g(to) on dil \ K(to), and the surface energy, for any fracture K D K(to) is 
proportional to Kri l {K), where T-L 1 denotes the one dimensional HausdorfF measure 
and k is a constant which is known as the toughness of the material. Here the 
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matrix A which appears in the integral in (1) is (fj,/2)Id, however in the paper we 
will also address the case of more general matrices A(x), which will be assumed to 
be uniformly elliptic and spatially Holder-continuous. 

The proof of existence for a crack K(t) satisfying the propagation criterions of 
brittle fracture as postulated by Francfort and Marigo in [10], was first proved 
by Dal Maso and Toader [8] in the simple 2D linearized anti-plane setting, then 
extended in various directions by other authors [4, 7, 11, 1]. 

In this paper we will freeze the "time" at a certain fixed value to, ancl therefore 
do not really matter exactly which model of existence we use. We will only need to 
know that such fractures exist, as a main motivation for our results. 

In the quasistatic model, the fracture K(t) is in equilibrium at any time, which 
means that the total energy cannot be improved at time to by extending the crack. 
Precisely, for any closed set K D K(to) such that K (to) U K is connected, and for 
any u g H 1 (Q\ (K(to) U K)) satisfying u — g (t ) on dQ \ (K(to) U K), one must 
have that 

E(t ) + nU l {K{t )) < I (AVu) ■ Vw dx + ^(K). 

Jn\(K(t )UK) 

This implies that the propagation of the crack is totally dependent on the external 
force g, and a necessary condition for K to propagate is that of the first order limit 
of the bulk energy, namely 

f0 , v E(to + h)- E(tp) 

(2) hmsup , 

to be greater or equal to k. The limit in (2) can be interpreted as an energy release 
rate along the growing crack, which is the central object of many recent works 
[6, 5, 3, 13, 14]. 

In all the aforementioned papers, a strong regularity assumption is made on the 
fracture K(t): it is assumed to be a segment near the tip in [5, 3, 13]; to our 
knowledge the weakest assumption is the C 1 ' 1 regularity in [14]. The main reason 
for this is the precise knowledge of the asymptotic development of the displacement 
u near the tip of the crack, when it is straight. Indeed the standard elliptic theory 
in polygonal domains (see e.g. Grisvard [12]) says that in a small ball B(0,e) (we 
assume the crack tip is the origin), if u denotes the minimizer for the problem (1), 
then there exists u € H 2 (B(0, e) \ K(to)) such that 

(3) u = Cy/rs,m(6/2)+u, 

(in polar coordinates, assuming the crack is {8 = ±7r}). In fracture theory, the 
constant C in front of the sinus is usually referred as the stress intensity factor 
(SIF). In [14], G. Lazzaroni and R. Toader proved that (3) is still true if K(to) is a 
C 1 ' 1 regular curve, up to a change of coordinates, and they base their study of the 
energy release rate upon this fact. 

The main goal of this paper is to extend (3) to fractures that are merely closed 
and connected sets, and asymptotic to a half-line at small scales. (We will need 
the technical assumption that the Hausdorff density is 1/2 at the origin, that is, 
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the length in small balls is roughly the radius — which basically means that K(to) 
admits a tangent, up to suitable rotations.) Our main result is as follows: 

Theorem 1.1. Assume that K := K(t ) C £1 C M. 2 is closed and connected, and let 
u be a solution for the minimizing problem in (1) with some a-holderian coefficients 
A : ft — >• S 2x2 . Assume that xq € K PI fi is a point of density 1/2, that is, 

H l (KC\B(x ,r)) 1 
lim sup i = - 

r^o 2r 2 

and that A(xq) = Id. Then the limit 

(4) lim - / (AWu) ■ Vu dx 

r J B (xo,r)\K 

exists and is finite. Moreover denoting Co the value of this limit, considering R r a 
suitable family of rotations, and taking 



/ 1C 

g(r, 6) := u(0) + \j sin(0/2), (r, 9) e [0, 1] x [-tt, tt], 

then the blow up sequence u r := r~^u(rR r (x — Xq)) converges to g and Vu r con- 
verges to both strongly in L 2 (B(0, 1)) when r — > 0. 

If A(xq) Id we obtain a similar statement by applying the change of variable 
x h-> y / A(xo)x (see Theorem 4.2). We also stress that a rigourous sense to the value 
of it(0) has to be given, and this will be done in Lemma 4.1. Besides, the exact 
definition of the rotations Rr will be given in Remark 2. 

Theorem 1.1 is a first step toward understanding the energy release rate for non- 
smooth fractures, and study qualitative properties of the crack path. It provides 
also the existence of a generalized stress intensity factor, that we can define as 
being the limit in (4), and which always exists without any regularity assumptions 
on K(to) of that of being closed and connected (see Proposition 7). 

Our main motivation is the study of brittle fracture, but of course Theorem 1.1 
contains a general result about the regularity of solutions for a Neumann Problem 
in rough domains, that could be interesting for other purpose. 

The proof of Theorem 1.1 will be done in two main steps, presented in Section 3 
and Section 4, which will come just after some preliminaries (Section 2). The first 
step is to prove the existence of limit in (4). For this we will use the monotonicity 
argument of Bonnet [2], which was used to prove existence of blow up limits for 
the minimizers of the Mumford-Shah functional. We will adapt here the argument 
to more general energies as the one with coefficient A(x), and also with a second 
member /. Notice that when / = we need only K to be closed and connected, 
whereas when / ^Owe need furthermore that K is of finite length. 

The second step is to prove the convergence strongly in L 2 of the blow-up limit 
u r := r~^u{rR r (x — xq)) and its gradient, to the function y^r sin(#/2). This is the 
purpose of Theorem 4.2, and the existence of limit in (4) is the first step, because it 
implies that Vu r is bounded in L 2 (B(0, 1)) which helps us to extract subsequences. 
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Notice that Bonnet [2] already had a kind of blow-up convergence for u ri ana- 
logue to ours in his paper on regularity for Mumford-Shah minimizers. The main 
difference with the result of Bonnet, is that here the set K is any given set whereas 
for Bonnet, K was a minimizer for the Mumford-Shah functional, which allowed 
him to modify it at his convenience to create competitors and prove some results on 
u. Here we cannot argue by the same way and this brings some interesting technical 
difficulties in the proof of convergence of u r . 



2. Preliminaries 

Let fi C R 2 , K C fl be a closed and connected set satisfying % 1 {K) < +oo 
(here V. 1 denotes the one dimensional Hausdorff measure), / <E L°°(f2), A > and 
g e ff^Q) n L°°(0). If K and K' are two closed sets of M 2 we will denote the 
Hausdorff distance by 



<Ih{K,K') := max I sup dist(a;, K'), sup dist(x,K) 

\x£K xeK' 

We also consider some a-H61dcr regular coefficients A : x A(x) e 6> 2x2 , 
uniformly bounded and uniformly coercive (with constant 7). We will use the 
following series of notations 

\\X\\ A :^XAX = (AX) ■ X = \\VAX\\ Id = \\yfAX\\. 

For simplicity we will assume without loss of generality that k = 1. We consider a 
slight more general energy than the one in (1) with a second member /, namely 

(5) F(u):= [ \\Vu\\ 2 A dx+\ [ \\u-f\ 2 . 

Jn\K A Jo. 

We will also allow the case A = and then we ask also / = and F is simply 



F(u) := [ \\ViL\\\dx. 



We consider a minimizer u for F among all functions v G i7 x (f2 \ K) such that 
v = g on d£l, i.e. u is a weak solution for the problem 

{Am - div^Vu = / in £1 \ K 
(AWu) ■ v = onK 
u = g on d£l 

It is well known that such a minimizer exists and is unique (up to additional constant 
if necessary in connected components of \ K when eventually / = 0), which 
provides a week solution for the problem (6). 

We begin with some elementary geometrical facts. 

Proposition 1. Let K C M 2 be a closed and connected set such that 
(7 hm sup = - . 
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For all r > small enough, let x r he any chosen point in K l~l dB(xo, r). Then we 
have that 

(8) lim -d H (Kr\B(xo,r), [x r ,x ]) = 0. 
r->o r 

Proof. Since K is closed, connected and not reduced to one point (because of (7)) 
we have that KndB(xo, r) is nonempty for all r small enough. Moreover since K is 
connected, there exists a simple connected curve T r G K that starts from xo and hits 
dB(xo, r) for the first time at some point y r € K ndB(xo, r). Since T r is connected 
we have that H 1 ^) > V}{[y r , x ]) = r and using (7) we get ft 1 (IV) < r + o(r). 
From the last two inequalities, since T r is a connected curve, it is then very classical 
using Pythagoras inequality to prove that 

(9) d H (T r ,[y r ,xo]) = o(r). 

Indeed, let z be the point in T r of maximal distance to [y r ,xo], and let h be this 
distance. Now let w be a point at distance h to [y r , xq], whose orthogonal projection 
onto [j/n^o] is exactly the middle of [y r ,xo]. Then the triangle (y r ,XQ,w) is isocel, 
and in particular minimizes the perimeter among all triangle of same basis and same 
height. Therefore, 

2y/(r/2) 2 + h 2 = \w - x r \ + \w- x \ < \z - x r \ + \z - x \ < H l (T r ) <r + o(r) 

which implies that h = o{r) and proves (9). 
Now (7) also implies that 

H 1 (KnB(x ,r)\T r ) = o(r), 

from which we deduce that 

sup{dist(x, T); x e K f] B(x , r)} = o(r) 

which implies dn{K n B(xo,r), [y r ,%o]) = °( r )- Finally (8) follows from the fact 
that dist(av, yr) = o(r) for any other point x r G K PI dB(xo, r). □ 

Remark 1. The density condition (7) does not imply the existence of tangent at 
the origin. One of such example can be found in [5, Remark 2.7.], as being a curve 
with oscillating tangent at the origin: exp(— t 2 )(cos(t)ei +cos(£)e2), t <E [0, +oo]. A 
further example is given by some infinite spirals, that turns infinitely many times 
around the origin but has finite length, and even density 1 /2 at the origin (thus is 
arbitrary close to a segment). To construct such an example one can consider the 
curve -f :t^ te w ^ where 6{t) G R satisfies lim^o d(t) = +oo and lim t ^ t0'(t) = 
(uniformly) like for instance 0(t) = sj\\\(— t). Then if K := 7([0, 1]) we have 

V}{K n B(0, r)) = f W(t)\dt= f ^l + t 2 6'(t) 2 dt^r + o{r) 
Jo Jo 

as desired. 

Remark 2 (Definition of R r ). As noticed in the preceding remark, the existence of 
tangent, i.e. the existence of a limit for the sequence of rescaled set -(K(~)B(xo, r)~ 
Xo), is not always guaranteed by the density condition. On the other hand if R r 
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denotes for each r > 0, the rotation that maps x r on the negative part of the first 
axis, then R r (^(K D B(xq, r) — xo)) converges to a segment. In the sequel, i? r will 
always refer to this rotation. 

Remark 3. There exists some connected sets such that j-K fl 5(0, r n ) converges 
to some radius in -8(0, 1) for some sequence r n — >• 0, and such that j-K n 5(0, t n ) 
converges to a diameter for another sequence t n — > 0. Such a set can be constructed 
as follows. Take a sequence q n —> such that 

(10) q n +l/q n — > 

The idea relies on the observation that thanks to (10), while looking at the scale of 
size q n , that is, in the ball 5(0, g„), all the piece of set contained in 5(0,q rl +i) is 
negligible in terms of Hausdorff distance. Therefore we can build two subsequences, 
one at the scales q 2n , and the other one at the scales q 2n+ \, that will not be seen 
by each other. 




B(0,q N ), N even 
K looks like a diameter 



B(0,q N ), N odd 
K looks like a radius 



Figure 1 . A crack tip with two different limits along different subsequences 

To do so, we consider the points A n := (0, §<7n) on the second axis of M 2 and we 
define K, as being IP x {0} union of all horizontal diameters of B(A2 n +i, <72n), that 
are connected to the first axis by their left extremities. In other words, denoting 
ei = (1, 0) and e2 = (0, 1) the two unit canonical vectors of M 2 , 

K:=MTe 1 u( (J (Re x + A 2n+1 ) n B(A 2n+1 , q 2n )) U (J ([0, -g n +i]e 2 + Z 2n ), 



where Z 2n is the left extremity point of the segment (Mei + A 2n +i) H B(A 2n+ i, q 2n ) 
(which is actually the horizontal diameter of B(A 2n +i, q 2 n)), sec Fig. 1. 
Then it is easy to see that, in the Hausdorff distance, 



1 

Kf)B(x,q 2n ) 

<?2ra 



lei n 5(0,1) 



and 



<?2n+l 



-{KnU{x,q 2n+1 ) 



"eifl 5(0,1) 



as desired. 
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Remark 4. Notice that a consequence of Theorem 1.1 for the example exhibited 
in Remark 3 is the following curious fact: even if n B(x,r) has no limit when 
r — > 0, the limit of i J" B / 0r ) l|Vu!| 2 as r — > exists thus has same value Co for 
any subsequences of r. Now, since K has density 1/2 along the odd subsequence 
r n = l2n+i, applying the proof of Theorem 1.1 for this subsequence we infer that the 
limit of the blow up sequence r n 1 ' 2 u(r n x) converges to ^2C'or/ir sin(0/2). But now 
regarding the limit in the even scales, r n = q2 n , as K is converging to a diameter, a 
similar proof as the one used to prove Theorem 1.1 would imply that the blow up 
sequence is converging to the solution of a Neumann problem in a domain which 
is a ball, cut into two pieces by a diameter. This implies Co = (because of 
the decomposition of u in spherical harmonics), so that actually returning to the 
odd subsequence, for which K is converging to a radius, we can conclude that 

— 1/2 

r n u(r n x) must converge to as well. 

It is well known that any closed and connected set K is arcwise connected, namely 
for any x ^ y in K one can find an injective Lipschitz curve inside K going from x 
to y (see e.g. [9, Proposition 30.14]). This allows us to talk about geodesic curve 
inside K, that connects x to y, which stands to be the curve with that property 
which support has minimal length. 

Definition 2.1. We say that K is locally-chord-arc at xo if there exists a constant 
C and a radius tq such that for every r < tq and for any couple of points y and z 
lying on K l~l 8B(xq, r) the geodesic curve inside K connecting y and z has length 
less than Cr. 

Proposition 2. Let K C M 2 be a closed and connected set satisfying the density 
condition 

(11) limsup ^-U\K n B(.t , r)) = \. 

r _>o zr I 

Then K is locally-chord-arc at Xq. 

Proof. The density condition (11) together with the fact that K is closed and 
connected guarantees that K is non reduced to one point, contains xq, and that 
dB(xo,r) n K is nonempty for r small enough. Let ?*o > be one of this radius 
small enough such that moreover 

(12) W}{K n B(x ,r)) < (1 + Vr < 3r . 

Let now y and z be two points in K n dB(xo,r) for any r < r$ and let V <Z K 
be the geodesic curve connecting y and z. Then T is injective (by definition since 
it is a geodesic) and in addition we claim that Y C i?(xo,3r). Indeed, otherwise 
there would be a point igT \ B(xq, 3r) which would imply 'H 1 (rn B(xq, 3r)) > 4r 
(because y and z are lying on 8B(xQ,r)) and this contradicts (12). But now that 
T C B(x ,3r), condition (12) again implies that H 1 ^) < U l (K n B(x , 3r)) < 4r 
which proves the proposition. □ 

In the sequel we will need to know that a minimizer of F is bounded. 
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Proposition 3. Let K be closed and connected, u be a minimizer for the functional 
F defined in (5) with f G and A > 0. Then 

NU < miii ^ 1 ^ maxdl/lloo, Iblloo)- 

Proof. It suffice to fix M := (min(l, A)) -1 maxdl/Uoo, ||si|oo) and notice that the 
function 

w := max(— M, min(it, M)) 

is a competitor for u and has less energy. By uniqueness of the minimizer we deduce 
that u = w. □ 



3. Bonnet's monotonicity Lemma and variants 
In this section we prove the existence of the limit 

lim - I \\Vu\\ A dx, 

r J yfA(ff)B{x,r) 

for any x G when u is a minimizer of F . Of course when x G VL \ K this is clear 
by the interior regularity of solution for the Problem 6, and the value of the limit 
in this case is zero. Therefore it is enough to consider a point x G K. 

The case of harmonic functions is slightly simpler than the general case, and need 
no further assumptions on K than being just closed and connected. This direclty 
comes from [2] and [9] but we will recall the proof in Subsection 3.2, that follows 
a lot the approach of G. David [9, Section 47]. Then we will consider the case of a 
non zero second member / but still with the classical Laplace operator, and finally 
in a third section we will adapt all the proofs to more general second order operator 
of divergence form. 

We begin with some technical tools. 

3.1. Technical tools. We will need the following 2 versions of the Gauss-Green 
formula. 

Lemma 3.1 (Integration by parts, first version). Let K be closed and connected, 
u be a minimizer for the functional F defined in (5). Then for any x G and for 
a.e. r such that B(x,r) C Q it holds 

/ \\Vu\\ A dy= (f-Xu)udy+ u(AVu) ■ v dU l . 

JB{x,r)\K JB(x,r)\K JdB(x,r)\K 

Proof. If u is a minimizer, then comparing the energy of u with the one of u + tip 
and using a standard variational argument yields that 

(13) / (AVu)-V<pdx= I (Xu-f)(pdx 

must hold for any function (p G 1 (O \ K) compactly supported inside Q. Let us 
choose tp to be equal to ip e u(x), where ips{ x ) = 9e{\\ x \\) is radial, and g e is equal 
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to 1 on [0, (1 — e)r], equal to on [r + e, +oo[ and linear on [(1 — e)r, (1 + e)r\. 
Applying (13) with ip = ip e u gives 

(14) / (AVu) ■ uVipe + I {A\7u) ■ i/'eVu = I (Xu- f)ip e u. 
Jn\K Jn\K Jn 

It is clear that ip E converges to lB(x,r) strongly in L 2 (fl), which gives the desired 
convergence for the second term and last term in (14). Now for the first term, we no- 
tice that ip £ is Lipschitz and its derivative is equal a.c. to 2 e\\x\\ ^-B(x.r+e)\B(x,(i-E)r) 
so that 

/ (AWu) ■ uVi/j £ = — I (AWu) ■ u-^-r 

Jn\K 2e J (B{x,r+e)\B{x,{l-e)r))\K \\ x \\ 

which converges to J dB ^ x r ^ K (AVu)-u v dH 1 for a.e. r by Lebesgue's differentiation 
theorem applied to the L 1 function r h-> f gB , x r )\i<-(^^ u ) ' v dH 1 . □ 

The first part of the next Lemma comes from a topological argument in [9] (see 
page 299). 

Lemma 3.2 (Integration by parts, second version). Let K C f2 be closed and 
connected, x € K and tq > be such that B(x,ro) C f2. For all r G (0,7"o) we 
decompose dB{x,r) \ K = [jj^jM Ij(f) where Ij(r) are disjoints arcs of circles. 
Then for each j € J(r) there exists a connected component Uj (r) of Q \ (Ij (r) U K) 
such that 

dU j (r)\K = I j (r). 

Moreover if u is a minimizer for the functional F defined in (5), then for a.e. 
r £ (0,ro) and for every j € Ij{r) we have 

(15) / (AVu) -vdT-L 1 = [ (Xu - f)dx, 

where v is the inward normal vector in Uj, i.e. pointing inside Uj. 

Proof. By assumption K is closed, so that dB(0,r) \ K is a relatively open set in 
dB(0, r) which is one dimensional. Therefore we can decompose dB(Q, r) \ K as a 
union of arc of circles as in the statement of the Lemma, namely 

dB(0,r)\K= \Jlj. 

(we avoid the dependance in r to lighten the notations). Let us denote by U^~ the 
connected component of \ (K U Ij ) containing the points of B(0, r) \ K very close 
to Ij, and similarly UJ is the one containing the points of il \ (K U 5(0, r)) very 
close to Ij . Then there is one between Up , that we will denote by Uj , which satisfies 

(16) dUj\K = Ij. 

The proof of (16) relies on the connectedness of K (see [9] page 299 and 300 for 
details: in our case the connectedness of K implies the topological assumption 
denoted by (8) in [9] that is used to prove (16) (which is actually (14) in [9])). 
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Then we want to prove (15) by an argument similar to Lemma 3.1 applied in Uj. 
For this purpose we consider as before a radial function but we need to separate two 
cases: if Uj C B(0,r) then we take the same function tpe(x) = 9e(\\%\\) where g £ is 
equal to 1 on [0, (1 — e)r], equal to on [r + e, +oo[ and linear on [(1 — e)r, (1 + e)r]. 
Now if Uj C \ B(0, r) we define tp e := 1 - ip E . 

Then we take as a competitor for u the function tp = lg-.^ e , where Uj is the 
connected component of fl \ K containing Uj. Notice that this is an admissible 
choice, namely tp S JT 1 (n \ K) and tp = on dft. 

Applying (3.1) with tp = ljj.V'e gives 

(17) [ (AVu) ■ V^ e = [ (Xu - f)ip s . 

JUj JO, 

As in the proof of Lemma 3.1, it is clear that ip £ converges to li/j strongly in L 2 (fl), 
which gives the desired convergence for the right hand side term in (17). Now for 
the left hand side term, we use as before that ip £ is Lipschitz and its derivative is 
equal a.e. to ^efa-y ^-B(x,r+s)\B(x,(i-e)r) ( w hh the correct sign depending on which 
side of Ij lies Uj) so that 

/ (AVu) ■ VV> £ = ±f / (AVu) • 

Jn\K ze JUjn(B(x,r+e)\B(x,(l-E)r)) \\ x \\ 

which converges to Jj (AVu) ■ vdH 1 for a.e. r by Lebesgue's differentiation theorem 
applied to the L 1 function r H> f gB r x r \ n y .(AVu) ■ v dH . □ 

3.2. Monotonicity for harmonic functions. So we arrive now to the first mono- 
tonicity result. The following proposition is one of the key points in Bonnet's proof 
of the classification of global minimizers for the Mumford-Shah functional [2] (see 
also Section 47 of Guy David's book [9] for a more detailled proof with slightly 
weaker assumptions than [2]). The same argument was also used in Lemma 2.2. 
of [15] to prove a monotonicity result for the energy of a harmonic function in the 
complement of minimal cones in M 3 , but the rate of decay obtained by this method 
is sharp only in dimension 2. Notice also that a similar argument with the elastic 
energy (i.e. L 2 norm of the symmetric gradient) of a vectorial function u : — > M 2 
seems not to be working. Notice that in [9] the assumption l-l}(K) < oo is needed 
whereas K is not necessarily connected. Here we do not suppose ^(K) < oo but 
we ask K to be connected which is a stronger topological assumption but weaker 
regularity assumption than [9]. 

Proposition 4 (Monotonicity of Energy, the harmonic case). [2, 9] Let K be a 

closed and connected set and let u be a solution for the problem (6) with A = Id, 
f = and A = (therefore u is harmonic in Q\ K). For any point xq G K we 
denote 

E(r) := f \\Vu\\ 2 dx. 

Jb(x ,t)\K 

Then r t— > E(r)/r is an increasing function of r on (0,ro). As a consequence, the 
limit lim^o E(r)/r exists and is finite. 
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Proof. Let us rewrite here the proof contained in [9] and [2] since we want to gen- 
eralize it just after. We assume without loss of generality that xq is the origin. 
Firstly, it is easy to show that E admits a derivative a.e. and 

(18) E'(r) := f \\Vu\\ 2 dx. 

JdB(0,r)\K 

In addition E is absolutely continuous (see [9]). Therefore, to prove the monotonic- 
ity of v i — y E(r)/r, it is enough to prove the inequality 

(19) E(r) < rE'(r) for a.e. r < r , 

because this implies (E{r)/r)' > a.e. 

We will need Wirtingcr's inequality (see e.g. page 301 of [9]), i.e. for any arc of 
circle I r C dB(0,r) and for g £ W 1,2 (I r ) we have 



(20) 



I | 9 -m 9 |W<(^) 2 | \gfdH 1 



where m g is the average of g on I r and g' is the tangential derivative on the circle. 
The constant here is optimal, and is achieved for the unit circle by the function 
sin((9/2) on the arc of circle ] — tt,%[. 

The first Gauss-Green formula (i.e. Lemma 3.1) applied in 5(0, r) yields, for a.e. 
radius r, 

(21) / \\Wu\\ 2 dx= f updH 1 . 

JB(0,r)\K JdB(0,r)\K OU 

Now K is closed, so that dB(0, r) \ K is a relatively open set in dB(0, r) which is 
one dimensional. Therefore we can decompose <9-B(0, r) \ K as a union of arc of 
circles as in Lemma 3.2, namely 

dB(0,r)\K=\Jl j . 

Next, we apply Lemma 3.2 to obtain for each of those arcs Ij, 

(22) / |W=0. 



r dv 

Denoting by rrij the average of u on Ij we deduce that 



(23) / u^dH 1 ^ [ (u-m^dn 1 



Returning to (21) and plugging (23) we get 



(24) / \\^u\\ 2 dx < £ / \u- mj \ 

JB(0,r)\K j&J Jl J 



On 
dv 
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Then by use of Cauchy-Schwarz inequality and ab < i-a 2 + |6 2 we can write 





du 




/ \u-nij\ 


dv 


dH 1 < I 



< 




2e 

Using Wirtinger inequality and setting e = 2r we deduce that 

4r 2 f ( du\ 2 e f f du 





du 


/ \u-mj\ 


dv 





dH 1 



2e V dr J 2 . T \d 



< 



IV 



du\ 2 f ( du x 2 



,dr 1 Jt, \ dv 
(25) = r I ||Vu 



Finally summing over j, using (24) and (18) we get (19) and the proposition is 
proved. □ 

3.3. Monotonicity with a second member /. Now we start to prove some 
variants of Bonnet's monotonicity Lemma. If / is non zero, then we obtain a 
similar result but we need further assumptions on K to be of finite length and 
locally-chord-arc. 

Proposition 5 (Monotonicity, with a second member). Let u be a solution for the 
problem (6) with A = Id, A > and f,g€ L°° , and assume that K is a closed and 
connected set of finite length. For any point xq € K we denote 

E{r) := j \\\Ju\\ 2 dx, 

JB(x a ,r)\K 

and we denote by N(r) G [0,+oo] the number of points of KC\dB(Q,r). We assume 
in addition that K is locally-chord-arc at point xq. Then there exists a radius ro 
and a constant C depending only on ||/||oo ; llslloo o,nd the locally-chord-arc constant 
of K such that 

r i — ^ E^l + C P{r) 
r 

is an increasing function of r on (0,ro), where P{r) is a primitive of N(r). As a 
consequence, the limit lim r _>.o E(r)/r exists and is finite. 

Proof. The proof is similar to the one of Proposition 4. We want to prove that the 
second member / is just a perturbation under control which does not affect the 
limit of E(r)/r. Precisely, this time we will prove the inequality 

(26) E(r) < rE\r) + CN(r)r 2 for a.e. r < r , 

with N(r) G ^([0,^]). This implies that £(CP(r) + E(r)/r) > thus r i-> 
E{r)/r + CP{r) is increasing and this is enough to prove the Proposition because 
the limit of P(r) exists at 0. 
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We assume xq = 0. Observe that since K has a finite length, we know that 
§K(ldB(0, r) is finite for a.c. r £ (0, ro). Actually we will need to know a bit more. 
If N(r) denotes the number of points of K n B(0, r), by [9, Lemma 26. f.] we know 
that N is borel mesurable on (0,ro) and that 

(27) / N(s)ds<H 1 (KnB(0,t)). 



This will be needed later. For now, take a radius r a.e. in (0,ro) such that N(r) < 
+oo and decompose S r := dB(0, r) \ K into a finite number of arcs of circle denoted 
Ij, for j = l..N(r). Moreover since K is closed and connected, for each j there 
exists a geodesic curve Fj c K connecting the two endpoints of Ij. We denote Dj 
the domain delimited by Ij and Fj. Since K is locally-chord-arc at the origin we 
infer that \Dj\ < Cr 2 . Notice also that Dj corresponds to the set Uj of Lemma 3.2. 
The Gauss-Green formula (Lemma 3.1) applied in B(0, r) yields 

r r N{r) f i 

(28) / \\Vu\\ 2 dx= / (/- Xu)udx+ V / u^dH 1 , 

JB{0,r)\K JB(0,r)\K i=1 J Ij " v 

and applied in Dj (Lemma 3.2) gives 

f pdH 1 =± f (f- Xu)dx, 
Jlj ov J Dj 

the sign depending on the relative position of Dj with respect to dB(0, r). Denoting 
by rrij the average of u on Ij we deduce that 

(29) / u^-dU 1 = I (u-m i )^dH 1 ± [ m., (/ - An) dx. 
hi ov Jj. dv J D . 



Now since u is bounded it comes \rrij \ < C, and we also have J2j=i < CN(r)r 
Moreover / is also bounded thus returning to (28) and plugging (29) we get 



ou 



dv 



dn 1 . 



(30) / || Vu|j z dx < CN(r)r z + > / \u - m,j 

JB(0,r)\K J = 1 Jl, 

Then the same computations as for proving (25) (i.e. using Cauchy-Schwarz in- 
equality, ab < a 2 + 1 6 2 and Wirtinger) we obtain 





du 




/ \u-mj\ 


dv 


dV. 1 <r J 



|Vu|r<fa; , 

and after summing over j, (26) is proved, as claimed. □ 

3.4. The case of more a general operator. We consider now the general case 
with a-H61der regular coefficients A : x M- A(x) € £> 2x2 , uniformly bounded and 
7-coercive with 7 > 0. For any x € fi and r > we define the ellipsoid 



B A (x,r) := y/A(x)(B(x,r)). 
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Proposition 6 (Change of variable). Let u be a solution for Problem 6 in £1, that 
we assume to contain the origin, and fix Aq := A(0). Then uo 'Aq is the solution 
of a similar problem in (^/A~q) ~ 1 (Q,\K) with coefficient A := (\/ Aq)^ 1 o Ao (y/ Aq)^ 1 
instead of A. In particular A(0) = Id, and 

(31) / \\\7v\\\dx = { \\Vu\\ 2 A Act{^/A~ )- 1 dx 



B(0,r) A JB A (0,r)) 

Proof. Let u be a solution for Problem 6, and consider v := u a (\/Aq). Then since 
\/A~v is symmetric we have Vv(y) = (y/A~o~)Vu((y/A~o)(y)) and 

\Vu\\\dx = f ||\/AVii|| 2 da; 

\\VAVuo y/A^\\ 2 det(y/A^)dx 



||\/AV^ 1 Vv\\ 2 dct(y^)dx 



(32) = / \\Vv\\ 2 A det(y/A^)dx 

with A := {\J Aq)^ 1 o A o (^/Aq)^ 1 . Therefore if u is a minimizer for the functional 
F defined in (5), then v must be a minimizer for the functional 

(33) F(v) := / \\Vv\\ 2 A dx + ± f \Xv - f\ 2 , 



A 



with Cl := \/Ao 1 (£1 \ K) and / := / o ^/Aq. Finally, the same change of variable as 
the one used for (32) proves (31), which completes the proof of the proposition. □ 

Here is now the analogue of Proposition 4. 

Proposition 7 (Monotonicity of energy for general coefficients). Assume that K 
is a closed and connected set. Let u be a solution for the problem (6) with some 
a-Holder regular coefficients A, and with A = and f = 0. For any point xo G K 
we denote 



E{r) := / ||V<u||iefcc. 

Jb a {x .t)\K 

Then the function 

r^m {1 + C r^ 
r 

is nondecreasing. As a consequence, the limit lim r _ > o(^'( r ')/ r ) exists and is finite. 

Proof. We will use a third time a variation of Bonnet's monotonicity Lemma, i.e. we 
will follow again the proof of Proposition 4. Let us assume without loss of generality 
that xo is the origin. First of all, up to the change of coordinates x M> y/ A(0)x and 
thank to Proposition 6 we can assume without loss of generality that A(0) = Id. 
In this case BA(x,r) = P>id{x, r) = B{x,r). 

The Gauss-Green formula (Lemma 3.1) applied in B(0,r) yields 



(34) 



/ \\Vu\\ 2 A dx = y2 [ u(AVu) ■ v dH 1 , 

JB(0,r)\K ■ Jlj 
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where dB(0, r) \ K = Ujlj. On the other hand Lemma 3.2 gives for each j, 

(AVu) ■ v dx = 0. 

f 

Denoting by mj the average of u on Ij we deduce that 

(35) / m(AVu) ■ v dH 1 = [ {u-m^iAVuj-vdl-L 1 . 



1 

< — 

4r 



Thus 

(36) f \\Vu\\ 2 A dx <J2 [ \u-mj\\(AVu) ■ v\ dU 1 . 

JB(0,r)\K j=lA 

Then by use of Cauchy-Schwarz inequality, ab < -^a 2 + rb 2 , and Wirtinger we can 
write 

\u — rtii 

\u- m \ 2 +r \(AVu) ■ u\ 2 
Ji 3 

(37) < r[ \\7u- T \ 2 +r[ \{AX/u)-v\ 2 . 

Jij J i j 

Now we want to recover the full norm ||Vu||a from the partial norms |Vu • r| and 
|(AVu) • v\. For this purpose we write 

|Vu • t\ 2 = \{AWu) ■ t\ 2 + \(Id- A)Vu ■ t\ 2 + [2(AVu) • r][((Id - A)Vu) • r], 

and we notice that, by Holder regularity of A and 7-coerciveness we have (the 
constant C can vary from line to line) 

I (Id- A)Vu-t\ 2 < \\{Id- A)Vu|| 2 

< \\Id-A\\ 2 c{R2M2) \\Vu\\ 2 

< CV 2a ||Vu|| 2 

(38) < 7 Cr 2a \\Vu\\ 2 A , 
and 

2\AX7u ■ r{Id - A)Vu ■ t\ = 2\AVu-r\\{Id- A)Vu-t\ 

< 2||AVu||Cr Q ||Vum 

< Cr^miUIIVwIli. 

Therefore summing over j and putting all the estimates together we have proved 
that for r small enough, 

/ l|Vu|fcda; < r[ {\A\7u ■ t\ 2 + Cr a \\Vu\\ 2 A ) + r [ \(A\Ju)-v\ 2 

JB(0,r)\K JdB(0,r) JdB(0,r) 

(39) = r f \\AVu\\ 2 + Cr l+a f \\Vu\\ 2 A . 

JdB(0,r) JdB(0,r) 
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By Holder regularity of A we infer that 

I|Va||^ r2) <i + cW 2 , 

which implies 



\\AVu\\ 2 < / \\VA\\ 2 \\VAVv" 2 

dB(0.r) JdB(0,r) 



< (l + Cr Q / 2 ) f ||VAVu|| 2 = (l + Cr a/2 ) f \\Vu\ 

JdB(0,r) JdB(0,r) 



2 

A- 



3(0, r) JdB(0,r) 

Therefore, since E'(r) = J dB ^ r s \\Vu\W we have proved for r small enough, 

E{r) < (r + Cr 1+a ' 2 )E'(r), 

and we conclude with Lemma 3.3 below, applied with the exponent a/2 S (0, 1). □ 

Lemma 3.3 (Gronwall type, version 1). Assume that E(r) admits a derivative a.e. 
on [0,ro], is absolutely continuous, and satisfies the following inequality for some 
a e (0,1) 

(40) E{r)< (r + Cr 1+a )E'(r), Vre[0,r ]. 
Then the function 

r 

is nondecreasing. As a consequence, the limit \im r ^o(E(r) /r) exists and is finite. 
Proof. Observe that a primitive of l/(r + Cr 1+a ) is 

(41) / "~~ t~~ dr = In ( -) =: h(r). 

J r + Cr l + a \{Cr a + l)iJ 

Hence (40) yields that 

(E{r)e~ h ^y = (-h'(r)E(r) +E'(r))e- h(r > > 0, 

in other words, 

r 

is nondecreasing. Therefore the limit of E(r)(l + Cr a )~ jr exists when r goes to 
zero, and since (1 + Cr a )i converges to 1, we obtain the existence of limit also 
for E(r)/r. Now by monotonicity, this limit is necessarily finite since less than 
E ( r °> l\ + CVq ) ° which is finite for some rg fixed. □ 



ro 

We also have an analogue of Proposition 5 in the context of general coefficients 
which is the proposition below. 

Proposition 8 (Energy estimate for general coefficients and second member). Let 

u be a solution for the problem (6) with a-Holderian coefficients A, A > and 
f, g G L°° . and assume that K is a closed and connected set of finite length. For 
any point xq G K we denote 

E{r) := { \\Vu\\\dx. 

JB A (x Q ,r)\K 
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We assume in addition that K is locally- chord- arc at point xq. Then the limit 
lim r _j.o(£'( r )/ r ') exists and is finite. 

Proof. We follow the proof of Proposition 5, with the changes already used in the 
proof of Proposition 7. The main difference with the preceding propositions is that 
we arrive now to the inequality 

E(r) < (r + Cr 1+a/2 )E'(r) + CN(r)r 2 . 

then we conclude with the Lemma 3.4 below. □ 

Lemma 3.4 (Gronwall type, version 2). Assume that E(r) admits a derivative a.e. 
on [0, ro], is absolutely continuous, and satisfies the following inequality for some 
a € (0,1) 

(42) E{r) < (r + Cr 1+a )E'(r) + CN{r)r 2 , Vr e [0, r ], 

with N integrable on (0,ro). Then the limit lim r ^Q E{r)/r exists and is finite. 
Proof. Let us first find a particular solution of the inhomogcncous equation 

(43) G{r) = (r + Cr 1+a )G'{r) + CN(r)r 2 . 

For this purpose, recall (see Lemma 3.3) that the solutions of the homogeneous first 
order linear equation 

r + G7' i+Q 

are given by 

f( r ) = A- — r— j — , A 6 E. 

Then from the method of "variation of the constant" we deduce that a particular 
solution for equation (43) is 

G(r) = A(r) 
because is). In particular we have 



(Cr a + l) 1 /" ' 

with X(r) = -C J r N{t)(Ct a + 1)^ dt (notice that N(t)(Ct a + 1)^ is integrable 



(44) lim = 0. 

r->0 r 

Now let us return to E(r), which is assumed to satisfy (42). If we subtract G(r) in 
the equation (42) we get 

H(r) < (r + Cr 1+a )H'(r) , 

where H(r) = E(r) — G(r). Therefore we can apply Lemma 3.3 to H which gives 
the existence of the limit 

,• {H(r)\ 
hm — — < +oo, 



r->0 

and we conclude using (44). □ 
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4. Blow up 

Here we prove the second part of Theorem 1.1 concerning the blow up sequence. 
Before going on with blow up limits at the origin, we start with a rigorous definition 
of u(0). Indeed, let u be a solution for the problem (6) with g € L°°, (A = / = 0) 
or (/ e L°° and A > 0). We suppose that K is a closed and connected set satisfying 
the density condition (7) at 0. Let R r the family of rotations given by remark 2 so 
that r~ 1 R r {K n B(Q, r)) converges to the segment [— 1, 0] x {0} when r goes to 0. 
For any r small enough we define A r := R r ~ 1 (B((r/2,0) 7 r/4)) and 



w\ L u{x)dx - 



Lemma 4.1 (Definition of it(0)). The sequence m r converges to some finite number 
that we will denote by it(0). 

Proof. We begin with a discrete sequence r n := 2~ n ro for some ro small, n £ N. In 
particular we assume ro small enough to have 

(45) — f ||Vu|| 2 tf.T < C VneN, 

r n JB(p,r n ) 

for some constant C that surely exists thank to Section 3. Since r^ 1 R r (KnB(0, r)) 
converges to the segment [—1,0] x {0}, we are sure that for r small enough and 
for every n, the ball B n := R~^(B(r n /2,0),3r n /8) does not meet K and contains 
both A Tn and A rn+1 . We denote by m n the average of u on B„ . Applying Poincare 
inequality in B n yields 



1 



\ArJ 



(u — m n )dx 



^ 77— T / l"-"»nl < °— I 
\ArJ JB n r„ 



and the same for m rn+1 so that at the end 



(46) |m r „ -m r „ +1 | < C— [ \\Vu\\dx < C ( f \\Vu\\ 2 dx) <Crl/ 2 

r « Jb„ \jb. / 



because of (45). In particular this implies that m r „ is a Cauchy sequence, thus 
converges to some limit £ € R. Now if is any other sequence converging to zero, 
we claim that the limit of m rk is still equal to £. To see this it suffice to find a 
subsequence r„ fc of r n such that r nk /2 < r^ < r nk and compare m rk with Tir„ fc by 
the same way as we obtained (46) and conclude that they must have same limit. □ 

Remark 5. In the future it will be convenient to introduce another type of averages 
on circles, namely 

■i- / udH 1 , 

A r Ja t 



m r 



with 



A r :=B((r,0), J)naB(0,r) 



It is easily checked that the sequence of rh r are aslo converging to u(0), i.e. has 
same limit as m r . 



SIF FOR NON-SMOOTH FRACTURES 



19 



We are now ready to prove the last part of Theorem 1.1. 

Theorem 4.2 (Convergence of the blow-up sequence). Let u be a solution for the 
problem (6) with g G L°° , (A = / = 0) or (/ g L°° and A > 0). We suppose that K 
is a closed and connected set satisfying the density condition (7) at the origin. We 
denote by u(0) the real number given by Lemma (4.1). Let R r the family of rotations 
given by Remark 2 so that r~ 1 R r (KnB(0, r)) converges to some segment £o when r 
goes to 0. If (r, 9) are the polar coordinates such that (t^A^O)) -1 (So) = (R _ x {tt}) 
we denote by Vq the function defined in polar coordinates by 



v (r, 6) := ^^ s in(0/2). 

Then 

u r := r-^{u{rR- l x) - u(0)) ^ v a o ^A(O), 
where the constant C$ is given by 



Co = lim ( \ / \\Vu\\ 2 A dx\ , 

and i/ie convergence holds strongly in L 2 (B(0, 1)) /or 6o£/i u r and Vu r . 

Proo/. We know that K r := ±R r (K) converges to the half-line R x {0} locally 
in R 2 for the Hausdorff distance. To simplify the notations and without loss of 
generality, in the sequel we will identify u with no i?" 1 , and K with R r (K) so that 
we can assume that i? r = Id for all r. We can also assume that u(0) = and as 
before, it is enough to consider the case when A(0) = Id because the general case 
follows using the change of variable of Proposition 6. 

As in the proof of Lemma 4.1, for any r we denote by m r the average of u on 
the ball B((r/2, 0), r/4). Then we consider the function u r (x) := r~*{u{rx) — m r ) 
defined in ±(Q\K). The domain l(Cl\K) converges to R 2 \K with K := R~ x {0}. 

We will prove that u r converges, in some sense that will be given later, to function 
in R 2 \ Kq that satisfies a certain Neumann problem. In the sequel we will work 
up to subsequences, but this will not be restrictive in the end by uniqueness of the 
limit. 

The starting point is that Vu r is uniformly bounded in L 2 (B(Q, 2)) (we start 
working in B(0, 2) for security but the real interesting ball will be B(0, 1)). Indeed, 
V« r (a;) = \frVu{rx), 

\\X7u r \\ 2 dx = I r\\\7u(rx)\\ 2 dx = - [ \\ Vu{x)\\ 2 dx. 

B(0.2)\K r JB{0.2)\K r r J B(0,2r)\K 

From Proposition 4, we know that i Je(Or) ll^' u ( a:: )ll converges to Co and we 
deduce (using the coerciveness of A), that Vw r is uniformly bounded in L 2 (B(0, 2)). 

Therefore we can extract a subsequence such that Vw r converges to some h, 
weakly in L 2 (B(0,2)), and 

(47) / \\h\\ 2 < liminf / ||VM r || 2 (i.T < C. 

JB{0,1) Jb(o,i) 
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Next we want to prove that in compact sets of 5(0, 2) \ Kq, the convergence is much 
better. For this purpose we introduce for any a > 

U(a) := {x e B(0,2);d(x,K ) > a}. 

The sequence u r is uniformly bounded in H 1 (U(a)) for any a. Therefore taking a 
sequence a n — >• 0, extracting some subsequence of u r and using a diagonal argument 
we can find a subsequence of u r , not relabeled, that converges weakly in H 1 and 
strongly in L 2 in any of the domains U(a). In other words, this subsequence u r 
converges weakly in 5^(5(0,2) \ Kq) and strongly in L 2 OC (B(0,2) 7 \K ) to some 
function uq <G -ff/ oc (i?(0, 2) \ Kq). By uniqueness of the limit we must have that 
Wuq = h a.e. in B(0, 2) and therefore (47) reads 

(48) / ||Vw || 2 < liminf / \\Wu r \\ 2 dx < C. 
Jb(o,i) Jb(o,i) 

Now wc want to prove that uo is a minimizer for the Dirichlet energy, and at 
the same time prove that the convergence hold strongly in L 2 (B(0, 1)) both for u r 
and Vu r . To do this we consider any function v £ H 1 (B(Q, 1) \ Kq) with v = uq in 
B(Q, 1) \ B(0, 1 - 5) and v = in 5(0, rf), for some small 6 > 0. The family of all 
such functions v is dense in the space of functions of H 1 (B(Q, 1) \ Kq) with trace 
equal to uo on dB(0, 1) \ Kq and therefore to prove that uq is a minimizer, it is 
enough to prove that 

f ||Vu || 2 <fa < f \\Vv\\ 2 dx 
Jb(o,\) Jb{q,i) 

for all such functions v. 

We denote by N r (s) the number of points of K r n dB(Q, s). As already used 
before, since by assumption H 1 (isr r H B(0, 1)) converges to 1 and 

1< [ 7V r (s)ds<H 1 (A" r nB(0, 1)), 
Jo 

we can extract a subsequence such that N r (s) — > 1 a.e. Then Fatou's lemma yields 

(49) [ liminf/" WVu^dH 1 ds < \immi [ [ WViirfdH 1 ds = G\, 

JO r JdB s r Jo JdB„ 

where C\ is closely related to Cq. This will allows us later to find a good radius s 
for which both N(s) ~ 1 and J gB \\ Vu^dH 1 is uniformly bounded. 

At this stage we only know that Vu r converges weakly in L 2 to Vuo. On the 
other hand, up to a further subsequence, we can find a measure [i such that |VM r | 2 (ia; 
weakly-* converges to \i. Let x € B(0, 2), p > such that B{x, p) C 5(0, 2) \ Kq. 
Let i/j be a smooth cutoff, with support in B(x,p), and equal to 1 in B(x,p/2). 
Then wc can write that 

(50) / (A r Vu r ) ■ X7(ip(u r — uq)) + r 2 Xu r (u r — uo)ip — r 3 ^ 2 f r (u r — uq)i/j = 
Jb{x, p ) 

where A r (x) = A(rx), f r (x) = f{rx) — Xm r , and (taking the limit in the "first" 
u r while freezing the test function (u r — Uq)iJj, and using the weak convergence in 
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H 1 {B(x 1 p) of u r to uq): 

(51) / (Vuo) • V(V>K - uo)) = 0. 

JB(x,p) 

Taking the difference of (50) and (51), and using the fact that u r —> uq strongly in 
L 2 (B r ), Vu r is uniformly bounded in L 2 (B r ) 2 , and A r — > Id uniformly, we obtain 
that 

lim / ||Vu r - Vu \\ 2 dx = 

r ~> a Jb(x,p/2) 

so that clearly, (B(0, 2) \ K ) = |V«o|| 2 dx: if \i has a singular part it must be 
concentrated on K . Moreover, we have s,0)}) = for all s G [0,2) but a 

countable number. (Observe that using any other test function in (50) and passing 
to the limit, we easily deduce that uq is harmonic in B(0, 2) \ Kq, but this will 
also be a consequence of the minimality of the Dirichlet energy which will soon be 
shown). 

Now from (49) we may choose s, 1 — S < s < 1, so that //({— s, 0}) = 0, N r (s) = 1 
for all r large enough, and 



liminf / WVurfdH 1 < 
r Job, 



>dB s 

In particular, upon extracting a further subsequence, we may assume that 

sup / ||Vit r +oo. 

r JdB s 

Then, by Sobolev's embedding, and using the fact that the averages fh r are uni- 
formly bounded (see Remark (5)), we deduce that there exists C > such that 

ll u r'||l,"=(aB s ) < C- 

We now consider any constant M > C and define 



A I 



(-M V (u r A M)) 



we have that uf, 1 — > u^ 1 in Lf oc (B(Q, 1) \ Kq), where u^ 1 is naturally defined as 
being Uq 1 := {—My (uq A M)). Up to a subsequence the convergence holds almost 
everywhere. But now, since the functions arc uniformly bounded, it converges also 
strongly in L 2 {B(0, 1) \ K ). 

Now, from the original function v G H 1 (B(0, 1) \ Kq), we want to construct a 
function v r e H 1 (B(0, 1) \ K r ) not much different from v. We denote by C^r the 
connected components of (B(0, 1) \ K r ) n {x < 0} containing (—1/2, ±1/2) and we 
define v r (x,y) as follows. In B(0, 1) n {x > 0} we set v r (x,y) = v(x,y). 



v(x,y) if y > 
w(x, — y) otherwise. 

if y < o 

u(x, —y) otherwise. 

And finally v r = everywhere else (i.e. in B(0, 1) n {a; < 0} \ (C+ n C~)). Then 
it is easy to see that v r G H 1 (B(0, 1) \ Jf r ), converges strongly to u in L 2 and 



InC+, v r {x,y) 
lnC~, v r (x,y) 
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t-B(p,i)\K r ^ v r converges strongly to 1b(q,i)\k Vv hi L 2 (B(0, 1)). However, by this 
procedure the trace on dB(0, 1) is not necessarily preserved. 

To get rid of that we let e < s — (1 — 5) , we pick a smooth cut-off ip s with compact 
support in B s , < < 1 and i/j e = 1 in i? s _ e , and we let 

v; = ip s v r + (1 - i> e )uf 

which converges strongly in L 2 to ip £ v + (1 — ip e )uQ as r — > 0. Next we write, since 
vf. = ixf = u r on dB s and u r is a minimizer, 

{A r \7u r )-\7u r +Xr 2 uj~2r 3/2 f r u r dx < [ (A r \7v e r )-\7v e r +\r 2 {v e r ) 2 -2r 3/2 f r v E r dx . 



Recall that |/ r | < C and |w r | < Cj\fr (by definition) so that 2r 3 / 2 f r u r = o(r) 
and we also easily check that 

Xr 2 (v e r ) 2 -2r 3/2 f r v E r dx = o(r) 
hence we focus on the other terms: we write for 5 > small, 



/ (A r \7u r ) ■ Vit r dx 
Jb, 



'B s 

< (l + v) [ ip 2 (A r Vv r ) • Vv r dx +C'/r] [ \\Vi> e \\ 2 \v r -u™\ 2 dx 
Jb b J B s 

+ C'/T] [ (l-^) 2 \\Vu r \\ 2 dx + o(r) 
Jb s 

Then sending r — > we obtain 

IIVuqII 2 < (1 + rfi [ i> 2 \\Vv\\ 2 dx 
Jb 3 

+ C'/r) [ \\V1> e \\ 2 \v-itf\ 2 dx + C'/Tin(B e (-8,0)) 
Jb s 

but on the support oiWip, v — uff is equal to (1 — ipc)^^ — uq). Therefore letting 
M tend to +oo we get 



HV^oll 2 < (1 + 7?) / i> 2 \\Vvfdx + C'/ W (B e (-s,0)) 
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finaly letting e — > 0, then r\ — > 0, and adding the integral over B(0, 1) \ B(0, s) on 
both sides (where v and uq actually coincide) we get the desired inequality, namely 



/ ||Vu || 2 ^ < / \\\7v\\ 2 dx 

JB(0,1) JB(0,1) 



j(o,i) Jb(o,i) 
which proves that uq is a minimizer. Moreover taking the particular choice v = uo 
in the same argument as before would give 

limsup / ||Vu r || 2 da; < / IjVwoll 2 ^ 

r->0 Jb(0,1) JB(0,1) 

and this toghcthcr with (48), implies the convergence of norms, which by the weak 
convergence yields the strong convergence in L 2 for the gradients, as desired. 



SIF FOR NON-SMOOTH FRACTURES 



2:s 



Finally all that we did in 5(0,1) could be done in any 5(0, i?) for R as large 
as we want, which gives a definition of uo in M 2 \ Kq. Moreover uo is of constant 
normalized energy. In other words we claim that s n- i J B ^ Q ^ ||Vuo|| 2 is constant 
in s, identically equal to Cq. Indeed, by the strong convergence in I? of Vit r , the 
value of \ J B( « ||Vu || 2 is given by 



lim — 

i — >o s 



\Vu r \ 



'B(O.s) 

which we actually claim to be equal to Co : a change of variable gives 

1 

B(0,s) r JB(0,rs) 



|Vd| 



1 



((Id- A)Vu,Vu). 



(52) 

>B(0,rs) '' JB(0,rs) 

The first term in (52) converges to sCq and the second term converges to zero 
because less than — ^4||l°°(b(o,j)) times something bounded. 

The latter implies that uq is the cracktip function. More precisely, we claim now 
that 



(53) 



u 



2C r 



sin(0/2). 



We shall give two different arguments for (53). The first one is very nice and 
due to Bonnet: returning to the proof of the monotonicity Lemma applied to ito, 
which says that s i-> - J B , Q , ||Vito|| 2 must be increasing (Proposition 4), since 



B(0,s) 



'B(O.s) 

jVuo|| 2 is actually constant in s, all the inequalities in the proof are 
equalities. In particular u must be the optimal function in Wirtinger inequality, 
thus it is the famous Cy/r sin(0/2) function. 

The second argument is to decompose uq in spherical harmonics, i.e. as a sum 
of homogeneous harmonic functions in the complement of the half line Kq, which 
Neumann boundary conditions on K . Now using that s 4 j Jb(Os) H^ u o|| 2 is 
constant we can kill all the terms of degree different from 1/2 by taking blow-up 
and blow-in limits. This implies that uq must be homogeneous of degree 1/2, and 
from this information it is not difficult to deduce (53) by looking at uo on the unit 
circle. 

Then, the exact constant C := J in front of the sinus can be easily computed 
by hand with the formulas 



du 
dr 



1 du 
r~dd 



1 



2v^ 



and 



duo 
dr 



1 



It comes 



RCq = 



IVuoII 2 = 



'B(0,R) 





du 


i 




2 r R f 


f 

)B{Q,R) 


dr 


+ 


dr 


Jo J- 



r<2 

—drdO = C 2 R- 
4 2 
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Finally, originally u r was converging to 




sin(0/2) up to subsequences, but 



by uniqueness of the limit we conclude that the whole sequence converges to this 
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function and this achieves the proof. 
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